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Abstract 

In this MHD-model of the heliosphere, we assume a Parker-type flow, and a 
Parker-type spiral magnetic field, which is extrapolated further downstream from 
the termination shock to the heliopause. We raise the question whether the he- 
liopause nose region may be leaky with respect to fields and plasmas due to non- 
ideal plasma dynamics, implying a breakdown of the magnetic barrier. We analyse 
some simple scenarios to find reconnection rates and circumstances, under which 
the heliosphere can be an "open" or a "closed" magnetosphere. We do not pretend 
to offer a complete solution for the heliosphere, on the basis of nonideal MHD the- 
ory, but present a prescription to find such a solution on the basis of potential fields 
including the knowledge of neutral points. As an example we imitate the Parker 
spiral as a monopole with a superposition of homogeneous asymptotical boundary 
conditions. We use this toy model for x < —R where R = 100 AU is the distance of 
the termination shock to describe the situation in the nose region of the heliopause. 
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1 Introduction 



In the past several calculations concerning the role of magnetic reconne ction 
in the vicinity of the heliopause have been performed. Fahr et al . (1986) (and 
references therein) calculated reconnection probabilities and gave estimations 
for reconnection rates, but without complete reconnection solutions of the non- 
ideal MHD equations. Here, we present a plane model of the heliosphere in 
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Fig. 1. The heliospheric outflow in the vicinity of the stagnation point, similar to 
the classical Parker outflow. The scale is given in units of 100 AU 

the framework of stationary nonideal MHD. We assume that the tail direction 
is the direction of the (see Figure 1), which is also the direction of the 

asymptotical flow around the heliosphere. The direction of invariance is the 
^-direction. Results depend mainly on the tilt angle between the asymptotical 
magnetic field to the asymptotical plasma velocity. We concentrate on inves- 
tigating the region where the interstellar plasma encounters the solar wind 
plasma. 



2 Stationary nonideal MHD— flows in 2D 

2.1 Basic equations 

In our plane model of the heliosphere we make the following assumptions: 

• we take only stationary fields, i.e. d/dt = 0, and take: d/dz = 

• incompressibility as a substitute for the energy equation is assumed. In- 
compressibility here means that the density of a fluid element moving with 
the flow does not change in time meaning that the convective derivative 
^ := |£ + v ■ Vp vanishes. Using the mass continuity §f + V • (pv) = 
we thus can derive with ^ = — pV • v — which implies V • v = and 
v ■ Vp = 0. Then we introduce the auxiliary flow field or streaming vector 
w with w = ^fpv leading to V • w — 0. 

Therefore, the basic MHD equations are given by the following set 



V -B = 
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w: = ^/pv (2) 

V-w = (3) 

VxB = Hoj (4) 

-V0 e + E e z + —w x B = R (5) 

V Qw 2 +pj = -(VxB)xB-(Vxiij)xfiJ (6) 

where p is the thermal pressure and p + \w 2 = II is the Bernoulli pressure. 
<p e (x,y) is the part of the electric potential whose gradient gives the poloidal 
field components (i.e. the components in x- and ^-direction) of the electric 
field, so that E = —V(j) e (x,y) + E e z , and E is a constant component of 
the electric field in the invariant z-direction. This is to ensure that V x E = 
— = (stationarity condition) and to fulfill = 0. R is the nonidealness 
and the other symbols have thei r usu al meaning. Eq. (5 ) is th e generalized 
Ohm's law (see IVasvliunasI f)l972h and lliesse fc Schindlerl (|l988l )). For w and 



B we make the following ansatz: 



w = V(xe 2 + w z e z =w p + w z e z (7) 
B = Va x e z + B z e z = B p + B z e z (8) 

where the index p on the right hand side of these equations stands for the 
poloidal part of the magnetic and the velocity field, and (, a, B z , and w z 
are functions of x and y. With these assumptions we find for the generalized 
Ohm's law from Eq. (5) 



w x B = (VC x e z + w z e z ) x (Va x e z + B z e z ) (9) 
which results in 

- (VC, e z , Va)e z - B Z V( + w z Va = sfp (^Va + ^VC + R- E e z ) (10) 



The momentum balance equation Eq. (6) can be written as 



VIT = AaVa + —{VB Z x e z ) x (Va x e z ) B Z VB Z 

Ho Ho Ho 

+ACVC + (Vw z x e z ) x (V( x e z ) + w z Vw z 

= --AaVa - — (VB Z , e z , Va)e 2 - —B Z VB Z 
Ho Ho Ho 

+ACVC + ftw g , e z , V()e z + w z Vw z (11) 
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The problem here is to get rid of the z-component on the right hand side of 
equation (11), as the ^-dependence of the Bernoulli pressure II has to vanish 
if the problem should be restricted to a two-dimensional problem. This can 
be solved by setting 



= {Vw z , e z , VC)e 2 



/'o 



-{VB z ,e z ,Va)e z 



= VC x Vw z Va x VB Z 

= (VC xVa)^-l (Va x VC 



da po 
which in general then leads to 



dB z 



(12) 

(13) 
(14) 



dw T dB z 
+ 



da d( 







(15) 



This is the condition for the vanishing of the z-component in the Euler- 
equation (11). This partial differential equation (Eq. 15) can be solved by any 
regular function S = S((, a) which defines the shear components as 



dS_ 

d( 

dS_ 

da 



B y 



/'o 



(16) 
(17) 



as long as both we are dealing with a real two-dimensional problem, and there 
does not exist an implicit function of C and a, or that C is an explicit function 
of a, so that we need to have 



0(C,Q) 
d(x,y) 



^0 



;is) 



almost everyw here in the consid ere d domain. Other cases has be en discussed 
for example by Tsinganosl ( 1981 ) or Goedbloed fc Lifschit d ( 1997 ) without the 
constraint of incompressibility but with the constraint of a vanishing nonideal 
term R. 

It should be mentioned, that due to the mass continuity equation for the 
stationary and incompressible flow the density p is a function of the stream 
function C, since from the mass continuity equation v ■ Vp = it follows 



y/P 



Vp 



w 



Vp = (VC x e z + w z e z ) ■ Vp 



d (p, C) 
d(x,y) 







(19) 
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which implies that 

p = p(0 



(20) 



Considering the last steps we can rewrite Ohm's law (Eq. (10)) in the form 

§f4 = VP(*'-£„) 

d(x,y) 

-B z w p = y/p + w p - 
w z B p = ^[R a + B p 



9C, 
da 



(21) 



where w p = \w p \, B p = \B P \. The first equation represents the z-component, 
and the other two equations the (- and a-components (i.e. the poloidal com- 
ponents) of Ohm's law. The contravariant and normalized components of R, 
which is represented by 



R = R< (VC)° + R a (Va)° + R z (Vz)° , 
are given by 



(22) 



R a 



R< 



R z 



R ■ 


[(v,)° 
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(VC) 


0" 


e z ■ 


;(vc)° 
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(Va 


>° 










0" 


R ■ 


"(Va)° 


X 






e z ■ 


"(VC)° 
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'1 










0" 
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(VC)° 
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e z ■ 


"(VC)° 
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(23) 



(24) 



(25) 



The superscript denotes the unit vectors in this direction. 

Next we look at the Euler or momentum balance equation (Eq. (11)): 



Vn = — -AaVa - —B Z VB Z + A(VC + w z Vw z 
Ho no 



(26) 
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with 



VP = — AaVa + A(V( 
Po 



1 ^2 1 , 



B 2 

P = 11 + -^B z - -w z = v H + -(v 2 + v 2y 



(27) 



(28 



This delivers the whole set of equations of motion to be solved 



OP 1 . 
-— = — Aq 



as 

#C 



«x 



— = -B M 
da Ho z 



(29) 
(30) 



with two potentials, which are explicit functions of the stream function ( and 
the (magnetic) flux funct ion a. These equations have already been derived by 
Neukirch fc Priest ( 1996f ) for the pure 2D case with B z = and w z = 0. 



2.2 Conditions for field line conservation 



We know from Vasvriunasl ( 1972h and Hesse fc Schindler ( 1988h that field line 
connection breaks down with respect to the plasma velocity Vp = v if, in our 
stationary case the equation 



V x (vx B 



XB 



(31) 



for an unknown function A, where v and B are given, cannot be fulfilled. For 
solutions of non-ideal MHD we have to find solutions with an appropriate 
function A. To get reconnection solutions it is necessary to violate condition 
(31). Therefore, within the first step, in order to fulfill Eq. (31) we take Ohm's 
law (9) and insert this into (31). 



A5 = Vx 



— (VC, V«, e z ) e z - -^-B Z V( + -^-w z Va 



.Vp 

1 d((,a 



Vp 
i 



\Vp d ( x >y) 



x i V ( — B z ) x VC + V y— w z j x Va 



.VP d(x,y) 



x e. 



.Vp 

_dl 

y/p da 



' 1 dB g \-. 
— — — Va x VC 
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I dw z w z dp\ - - 

/ 1 d{C,a) \ ^ ( I OB 1 dw z w z dp\ -> -> 

(32) 

For the poloidal part and the shear part B z of the magnetic field it follows 
immediately that 

— !£4=A(«) (33) 
(dB z dw z w z dp\ 

{-d^ + ^-YpTc) A{a) = X{a)Bz (34) 

The above relations can be shown to be valid in the following way: From 



y/p d(x,y) 



x e z = (35) 



one derives 



and taking the curl of both sides of the above Eq. (36) one gets 

= Vx(AVa) VA x Va = =>- A = A(a) (37) 

and defining 

AVa :=: VA(a) = A'(a)Va (38) 
we thus get Eq. (33) and Eq. (34). 

For the condition of a frozen in magnetic flux it is therefore necessary that 
A = which implies that A(ct) = constant. Equation (34) can be expressed in 
the form of an equation similar to the telegrapher's equation for the unknown 
shear function S 



d 2 S d 2 S ^dS T dS 

-^ + ^- D Tr L s- a =° (39) 
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where 



D = 2^f and 1 =-"»-*r i (40) 



It is possible to rewrite equation (33) with the help of the chain rule 
d ( 1 d ((, a) \ fd((,a)\ 1 f da d da d \ f 1 d (£, a) \ 



d( \ yfp d (x, y) J \d(x,y)J \dy dx dx dy J \v/p d (x,y) 



which is equivalent to 



(41) 



1 d(C,a)Y 1 (B p -V)f 1 d(C,<*)\_(B p -V)/ 1 a(C,a) \_ 0(12) 



\Vpd(x,y)J \B p \ \y/pd(x,y)J \B p \ \^/pd(x,y) 

In the following we will use Eq. (41) to define line violation and therefore 
reconnection rates. 

It should also be mentioned that the z-component of Ohm's law (Eq. (21)) 
indicates, that the asymptotical flow component within the plane around the 
heliosphere, which may be inclined with respect to the asymptotical magnetic 
field components within the plane, determines the constant electric field in the 
invariant direction. Consequently, if we restrict to B z = and w z = and if the 
plasma is ideal everywhere and the flow is field-aligned at infinity, then the flow 
is field-aligned everywhere. This is due to the fact that the electric field in the 
^-direction is invariant. So if the electric field in the invariant direction vanishes 
at infinity, the Jacobian matrix in Eq. (21) vanishes. This then implies that the 
flow is in fact field-aligned eve rywhere. A similar d iscussion for axissymmetric 



equilibria can be found e.g. in IContopoulo s (1996). 



2. 3 Application to a plane in the vicinity of the equatorial plane of the sun 



We set S = w hich results in B z = and w z 
( Parkerl (|l96lh :'or our investigation 



and take the Parker outflow 



C = w ooy + Do arctan 



(43) 



The Parker flow field is defined by a potential field. In this article we assume 
the auxilliary flow field and the magnetic field to be potential fields, i.e. A( = 
and Aa = 0. The stream lines of the Parker outflow given by the potential £ in 
Eq.(43) can be seen in Fig. 1. In general, potential fields can be represented by 
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Laurent series. We define A = m + ia as the complex flux function with <p m , 
the scalar magnetic potential, as the real, and a as the imaginary part of A. 
Using the assumption of a finite number of neutral points u k and asymptotical 
boundary conditions, lim B 

x,y—>oc 



-Bqo, we can write: 



N 

A = B^u + C In u + C ' 

u=l 



,U~ 



n( \ dA 
B{u k ) = — 

du 



u=u k 



The roots of the polynomial B(uk) = are correlated to multipole moments: 



//- ' + ^u N - £ ^u N ^ = J] (« - u k ) = (44) 




Boo v=\ Boo 

N+l 

=>C = -B oo y £u k Ci, = (-!)"— E I II u k , (45) 

where 1 < v < N. The same procedure is valid for potential flows given by 
A( = 0. Therefore, we can draw the conclusion, that the topology determines 
the global geometry of a heliospheric model, if the magnetic neutral point, i.e. 
stagnation point distribution for such a potential magnetic field is known. The 
field lines, crossing a magnetic neutral point are called magnetic separatrices. 
In the case of the heliosphere this magnetopause encloses the inner region of 
the heliosphere, therefore called heliopause. But on the other hand, there is 
also a 'hydropause', which is marked by the separatrix of the plasma flow. 
This separatrix passes through the stagnation point and leads to the following 
questions: what is the 'heliopause'? Is it the 'hydropause', or is it the 'mag- 
netopause' ? If hydro- and magnetopause are not identical, does this include a 
breakdown of the magnetic connectivity with respect to the bulk plasma flow 
and therefore magnetic reconnection? How large is the reconnection rate or 
line violation rate in the case of potential fields with an incompressible flow? 
This includes the question, whether the heliosphere is open or closed with 
respect to the counterstreaming interstellar medium. 



3 Reconnection solutions 



We will analyze the magnetic reconnection processes due to the simplified 
scenario we have developed in the last subsections with the h elp of the cr i tereon 
for line and flux conservation. This relation is given, e.g by Vasvliunasl ( 19721 ) 
and references given therein. Our assumpt ion here is, th at the magnetic field 
looks like a Parker -spiral (see e.g. theory: Parker ( 1958i l and measurements: 
Smith et al.1 (|l997h ) by imitating the azimuthal component and spiral pattern 
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Fig. 2. The heliospheric Parker spiral like magnetic field, with 10° (left panel) and 
40° (right panel) tilting of the asymptotical magnetic field with respect to the 
asymptotical classical Parker (out)flow; the asymptotical Parker flow is aligned to 
the x-axis. Scale as in Fig. 1 

of the field lines using a complex monopole moment representation. Hereby, 
the complex monopole moment is given by C = C 0r + iC i, which results in 
the following representation of a: 

Q[A) = a = B x>00 y - B V:00 x + C 0r arctan f + C 0i In (^J x 2 + y 2 ^j (46) 



Solving for the neutral points gives the relation between the monopole moment 
on one hand and the coordinates of the neutral point and the asymptotical 
components of the magnetic field on the other hand 

Cor = —B Xi00 xn (—By )00 y]sr) Coi = —By^XN (— -Re.ooZ/jv) (47) 



How does this given potential field influence the strength and importance 
of the resistivity R z 7 With the help of the asymptotical boundary condi- 
tions we get for the asymptotical electric field: lmx ri3/ _ >00 (il 2 — E ) — —E = 

We calculate reconnection rates (=line violation rates), by using the derivative 
of Ohm's law yielding 

d( d (x, y) 



where the Jacobian should be regarded as a function of C and a. We use Va- 
syliunas' expression of line conserving flows (|Vasvliunasl ()l972h V which shows 
the condition for line freezing. Line freezing is fulfilled, if r vanishes every- 
where. If r vanishes only in the ideal region x, y —>■ oo, line freezing is violated 
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Fig. 3. The line violation rate, with 10° (left panel) and 40° (right panel) tilting of 
the asymptotical magnetic field with respect to the asymptotical classical Parker 
(out)flow; the asymptotical Parker flow is aligned to the x-axis. Scale as in Fig. 1 

and reconnection is taking place. In this case we assume, that stagnation and 
magnetic neutral point are identical to fix at least one point of these two 
pauses. With increasing angle between aJ^ and B^, the Parker spiral winding 
grows, and also the deviation of the shapes of magneto- and hydropause (see 
Figs. 1 and 2). Fig. 3 shows that the line violation rate is much stronger in the 
lower right region, where the angle 7 between the asymptotical magnetic field 
and the flow is larger. Here the magnetic field lines of the Parker-spiral have 
a higher curvature in the case of larger inclination (see right panel of Fig. 3). 
In Fig. 3 in front of the heliopause nose there is a bulge which shows that the 
reconnection rate here is larger than at other locations in the vicinity of the 
stagnation point. 



4 Conclusions 

We presented a method for calculating reconnection rates for a driven recon- 
nection process. Further investigations need to consider non-potential fields to 
account for spontaneous processes including Ohmic heating. As we have shown 
it would be highly valuable to obtain observational knowledge of the location 
of magnetic neutral points and stagnation points in the heliosphere, because 
on the basis of that knowledge a unique solution could be given according to 
our prescription. 
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